OUTP-12-15P 



The impact of heavy-quark loops on LHC dark matter searches 

Ulrich Haisch, 1 '^ Felix Kahlhoefer, 1 '^] and James Unwin 1, 2 [^] 

1 Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 
1 Keble Road, Oxford 0X1 3NP, United Kingdom 
2 Mathematical Institute, University of Oxford, 24-29 St Giles, Oxford, 0X1 3LB, UK 

(Dated: August 24, 2012) 

If only tree-level processes are included in the analysis, LHC monojet searches give weak con- 
straints on the dark matter-proton scattering cross section arising from the exchange of a new 
heavy scalar or pseudoscalar mediator with Yukawa-like couplings to quarks. In this letter we cal- 
culate the constraints on these interactions from the CMS 5.0 fb -1 and ATLAS 4.7 fb -1 searches 
for jets with missing energy including the effects of heavy-quark loops. We find that the inclusion 
of such contributions leads to a dramatic increase in the predicted cross section and therefore a 
significant improvement of the bounds from LHC searches. 
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I. INTRODUCTION 

It is a remarkable fact that one of the leading exper- 
iments for dark matter (DM) detection is the LHC. Al- 
though any DM particles produced at the LHC will es- 
cape from the detector unnoticed, we may observe large 
amounts of missing transverse energy if a single 

jet (j) is produced in association with a pair of DM par- 
ticles. The experimental search for such monojet events 
provides model-independent bounds on the interaction 
strength of DM with quarks and gluons, constraining the 
same parameters as direct detection experiments [H [2] ■ 
These searches place the leading (and in some cases only) 
limits on models of DM over certain regions of parameter 
space. 

If the mediator of the DM interaction is sufficiently 
heavy, it can be integrated out to obtain an effec- 
tive higher-dimensional operator that describes the low- 
energy interactions between DM and Standard Model 
(SM) states. For instance, the interactions of a fermionic 
DM particle yb with SM quarks q via a heavy scalar me- 
diator $ with quark couplings proportional to the quark 
mass m q can be described by the four-fcrmion operator 



(1) 



and the DM pair production proceeds through tree-level 
diagrams like those shown on the left in Fig. [T] Unfortu- 
nately, the suppression scale A s in ([!]) is very difficult to 
constrain with LHC monojet searches [3TTT0] since the ini- 
tial state contains (apart from gluons) only light quarks, 
which by assumption have small couplings to the scalar 
mediator <E>. 

In this letter we observe that the situation changes 
dramatically beyond tree level, since now loop graphs in- 
volving virtual top quarks give rise to a j + $ T signal. 
Examples of such Feynman diagrams are displayed on the 
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right in Fig. [T] We find that by including these loop con- 
tributions, the predicted cross section for monojet pro- 
duction increases by more than two orders of magnitude 
leading to bounds on A s which are about a factor of 3 
stronger than the constraints from the simple tree-level 
analysis with DM-quark operators. To our knowledge, 
this important observation has not been made before. 

This paper is structured as follows: in Sec. [Hj we 
present the various operators which we intend to exam- 
ine. In Sec. |III| we calculate the constraints arising from 
the latest LHC searches for jets with f! T , including the 
full top-quark mass dependence of the squared matrix el- 
ements for DM pair + jet production. Our calculation 
is performed at the leading order (LO) in QCD, but we 
will comment on the importance of higher-order effects 
as well as the applicability of the heavy top-quark ap- 
proximation. Finally, in Sec. |IV| we discuss the impact of 
the enhanced monojet limits in the context of relic den- 
sity constraints and recent results from direct detection 
experiments. 



II. EFFECTIVE OPERATORS 

Our focus will be on interactions that result from the 
exchange of a new heavy scalar or pseudoscalar state 
which connects SM quarks to DM. In the case that the 
scalar mediator is a SM singlet it can couple to quarks 
via mixing with the Higgs and the induced couplings will 
be proportional to the SM Yukawa couplings: 



Tft — 
V 



(2) 



where v ~ 246 GeV is the Higgs vacuum expectation 
value. Integrating out the mediator $, we obtain the 
effective interaction ([!]) and the scale A s is related to the 
fundamental couplings by A;? = vM\j{g q g^). If DM is a 
scalar <f> rather than a fermion ip, we obtain in complete 
analogy the effective operator 



(3) 



2 




The other type of portal interaction for which the effec- 
tive operator may naturally have Yukawa-like couplings 
involves heavy pseudoscalar mediators. The correspond- 
ing effective operator is [TT] 



i 



"'q _ K 7 R / 
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This case is of particular interest, since direct detection 
signals associated to O^ are spin-dependent and sup- 
pressed by powers of q/mpj <C 1, where q is the momen- 
tum transfer and toat is the mass of the target nucleus. 
Thus, mono jet searches provide the only manner to ob- 
tain constraints on A p for the foreseeable future. 



III. LIMITS FROM MONOJET SEARCHES 

In their most recent analysis with an integrated lumi- 
nosity of 5.0 ft)" 1 at yfs = 7TeV the CMS collabora- 
tion pQ found no apparent excesses in their searches for 
jets with p T providing the leading monojet bounds on 
DM. CMS considered events with $ T > 350 GeV, pro- 
vided there was a primary jet with transverse mo- 
mentum pr > 110 GeV and pseudorapidity \t)\ < 2.4. A 
secondary jet (j 2 ) with p T > 30 GeV was also permitted 
if the two jets are not back-to-back: |A0(ji, j 2 )| < 2.5. 
Events with high-p^ tertiary jets, electrons or muons 
were vetoed. This null result excludes new contributions 
to the production cross section in excess of 0.032 pb at 
95% confidence level (CL). The ATLAS search [2] em- 
ploys very similar cuts and finds a comparable bound on 
the cross section. We have checked explicitly that us- 
ing the ATLAS data instead of the CMS data does not 
modify our results within errors. 

To calculate the predicted monojet cross section for 
the operators introduced above we have implemented 
each of them in FeynArts [12j and performed the com- 
putations with FormCalc and LoopTools |13j . Fur- 
thermore, as an independent cross-check, we have ver- 
ified our findings with MCFM 14 , modifying the process 
p + p — >• H (A) + j ^ t + t~ + j, which is based on the 
analytical results of [T3] for the scalar Higgs case and 
[16] for the pseudoscalar Higgs case. Both computations 
utilise MSTW2008L0 parton distributions [T7j. In our cal- 
culations we do not consider the effects of parton show- 
ering and hadronisation or the contribution of additional 
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TABLE I. Bounds on A at 95% CL from the CMS 5.0fb _1 
search for jets with $ T including loop-level processes in the 
analysis. The quoted errors reflect the scale uncertainties. 
See text for details. 



jets. As discussed in [HHH], the first two simplifications 
are justified, because the primary jet has sufficiently high 
Pt, which renders the impact of non-hard QCD radiation 
small. We will comment on the effect of two-jet events 
below. 

Our results are summarised in Tab. [T] and displayed in 
Fig. [2] for the case where DM is either a Dirac fermion 
or a complex scalar. For Majorana fermion or real scalar 
DM, the predicted cross sections are larger by a factor 
of 2 and so the bounds on A are stronger by roughly 
12%. For the scalar (pseudoscalar) operators we find 
that including the loop-level processes in the calculation 
increases the predicted monojet cross sections by a factor 
of around 500 (900), the precise value depending on the 
DM mass and the choices of renormalisation (hr) and 
factorisation (/zp) scales. These numbers translate into 
an increase of the bounds on A by a factor of 2.8 (3.1). 
Note that the limit on A s which we obtain from the tree- 
level processes alone is stronger than the one reported 
in [2] by about 25%, because we include bottom quarks 
in the initial state, which we find to give the dominant 
contribution at tree level. 

Before examining the impact of additional QCD radia- 
tion, we first discuss whether the above results could have 
been obtained without performing an actual loop calcula- 
tion, but rather by employing the heavy top-quark mass 
limit. For the operator Of, the effect of heavy-quark 
loops can be described in this approximation in terms of 
the effective operator 
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FIG. 2. LHC monojet bounds on the scale A at tree-level (red) and loop-level (green) for the effective operators Of, Op , and 
Of (from left to right). The black dashed curves indicate the requirement for the correct relic density. The width of the bands 
reflect the scale uncertainties. See text for details. 



The operator induced by Of (O^) is obtained by re- 
placing with 0t0 (G^G 11 ^ with Gl v G a ^). Effec- 
tive interactions like ^ have been studied previously in 
the context of monojet searches [3] [HI E] . In fact, for 
m t — > oo bounds on these operators can be translated 
into limits on Of'* 1 and O^ by the simple identifications 
A s = A 9 /(3tt) 1 /3 a nd A p = A 9 /(2tt) 1 / 3 . 

The scales involved in j + $ T production (i.e. the pr 
and the DM mass) are, however, not necessarily small 
compared to the top-quark mass, which implies that the 
infinite mass limit employed to obtain ^ is not a good 
approximation [19j . Numerically, we find that applying 
the mt — > oo limit overestimates the monojet cross sec- 
tions by a factor of around 3 for small DM mass and that 
the quality of the approximation rapidly degrades with 
increasing px cut and DM mass, resulting in errors of 
up to a factor of 40. These numbers imply that the cor- 
responding limits on A would be too strong by a factor 
of 1.2 (1.9) in the best (worst) case. This clearly shows 
that one cannot use Of g (or analogous operators) to infer 
faithful bounds on the quark-DM contact operators. 

To assess the theoretical errors in our analysis, we have 
studied the scale ambiguities of our results. We find the 
scale fj, which determines a s (fi) dynamically, i.e. we de- 
fine /i = £pt = fJ>R = [If and evaluate it on an event-by- 
event basis. Following common practice, we have varied 
£ in the range [|,2]. We find that the predictions of 
the tree-level and loop-level monojet cross sections cal- 
culated in this way vary by around ±20% and ±50%, 
respectively. The resulting uncertainties on A stay be- 
low ±2 GeV and ±15 GeV. They are given in Tab. [I and 
reflected by the width of the coloured bands in Fig. 2} 

The pronounced scale ambiguity of the loop-level re- 
sult indicates that next-to-leading order (NLO) correc- 
tions might be large. Indeed, the -ftf-factor represent- 
ing the ratio between the NLO and LO cross sections of 
p + p — > H + j evaluated in the m t —> oo limit, varies 
between 1.2 and 1.8 depending on the kinematic region 
and choice of parton densities [2"0Tf2"^| . Using MCFM, we 
explicitly verified that the same range of if-factors also 
apply to the monojet signal, although the imposed px 
cut is significantly higher in this case than in the case of 
Higgs + jet production. 



Another important and related issue is the relevance 
of events with a second high-py jet- Such events are al- 
lowed in the CMS analysis, as long as the two jets are 
not back-to-back. To estimate the cross section for events 
with two jets and $ Tl we use the cross section for H + 2j 
implemented into MCFM in the limit m t -> oo [23], [24]. 
This limit is known to work better in the case of H + 2j 
than for H + j [35] . We find that the dominant contribu- 
tion to two-jet events arises from processes that resemble 
monojet events but have an additional gluon in the final 
state. Allowing a second jet with large px increases the 
total cross section by about a factor of 2, consistent with 
the more accurate simulations from CMS pQ. 

The inclusion of NLO corrections and of two-jet events 
is hence expected to strengthen the bounds on A, pos- 
sibly by as much as 25%. However, we expect this im- 
provement to be somewhat smaller in the full calculation 
with resolved top-quark loops (in particular, for high px 
cut and large DM mass), so we prefer not to include 
these effects in our results and give a conservative bound. 
Clearly, a more careful analysis of finite top-quark mass 
effects at 0{a 4 s ) (along the lines of the recent Higgs + 
jet study 26j) is required to make a more quantitative 
statement. This study is left for future work. 

Finally, we note that the values of A that we can ex- 
clude with the current data are relatively low compared 
to typical LHC energies and thus we anticipate that the 
new scalar or pseudoscalar mediator may be produced 
on-shell in pp collisions, unless the couplings g q and gu 
in |2| are large (see [HI M In this case the analy- 

sis becomes more model-dependent and monojet bounds 
may be modified [28]. With increasing integrated lumi- 
nosity, however, the bounds on A will improve, making 
the method advocated here more reliable. 



IV. DISCUSSION 

The suppression scale A which can be constrained 
through monojet searches also enters in the formula for 
the DM annihilation cross section into quarks and the 
DM-proton scattering cross section. Consequently, lim- 
its on A can also be inferred from measurements of the 
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relic density and DM direct detection experiments. Here 
we compare these results with the bounds derived in the 
previous section. 

The annihilation cross sections for the three operators 
introduced in Sec. [il] are given e.g. in [10]. Using the 
expansion (av) — a + 6b/x + 0(l/x 2 ) with x = tudm/T, 
the DM relic density 51dm after freeze-out is given by [25] 



DM' 



1.07 x 10 9 

GeV 



Mi 



pi V-9 



(.. 



(6) 



where is the number of relativistic degrees of freedom 
at the freeze-out temperature Tf. The coefficients a and 
b are functions of A and moM- If we assume that the 
interactions between the DM particle and SM quarks are 
the dominant ones at freeze-out, we can determine A as 
a function of todm by the requirement that the observed 
DM relic abundance J7dm^ 2 — 0.1109 [30] is reproduced. 

The values of A which give the correct relic density are 
indicated by the black dashed curves in Fig. [2] From the 
intersections of the monojet bounds and the relic density 
constraints, we obtain the following lower bounds on the 
DM mass: 



Of 

ot 



> 173 (174) GeV 

> 60 (85) GeV , 
m > 1.5 (1.6) GeV . 



(7) 



Here the values in brackets apply if DM is a Majorana 
fermion or a real scalar instead of a Dirac fermion or a 
complex scalar. Note that the relic abundance require- 
ments may differ in alternative frameworks, for instance 
asymmetric DM (see e.g. [10]). Also, additional anni- 
hilation channels, e.g. into leptons, gauge bosons or new 
hidden sector states, may avoid an overproduction of DM 
and therefore help reduce tension with experiments. 

For the scalar operators Of and Of, the DM direct 
detection cross section is given by 



> - 



f4 m l f 

7T A 6 



A 4 n 



(8) 



where m p is the proton mass, \i v is the DM-proton re- 
duced mass and / ~ 0.35 is the scalar form factor of the 
proton [31]. The scattering cross sections for Majorana 
or real scalar DM are larger by a factor of 4. Utilizing 
the formulas ([8]), the monojet bounds on A can be trans- 
lated into limits on a pi which in turn can be compared 
to the exclusion limits obtained from XENON100 [52] . 
XENON10 [33] and CDMS-II (3J [35] as well as to 
the best-fit regions obtained from DAMA [36] and Co- 
GeNT [37] [35]. 

Our final results are shown in Fig. [3] for the case where 
the DM particle is a Dirac fermion. All shown bounds 
and best-fit regions represent 95% CL ranges. For large 
values of m^, as indicated by the relic density constraints, 
direct detection experiments give stronger bounds than 
monojet searches. For ~ 10 GeV, the bounds become 
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FIG. 3. LHC monojet bounds on the DM-proton cross sec- 
tion for the operator Of at tree-level (red) and loop-level 
(green) compared to various results from DM direct detection 
experiments. While the tree-level monojet bound is too weak 
to constrain the parameter regions favoured by DAMA and 
CoGeNT, the loop-level bound clearly excludes these regions. 



comparable, while below this value the bounds from LHC 
searches are far superior. We observe that the inclusion 
of loop-level processes gives a pertinent improvement of 
the monojet bounds, in particular because it excludes 
the possibility that the CoGeNT excess or the DAMA 
modulation arise from the interactions of a heavy scalar 
mediator. 

For the pseudoscalar operators the DM direct detection 
cross section is spin-dependent and suppressed by q 4 /m 4 , 
so that no relevant bounds on A can be obtained from 
direct detection experiments. Consequently, the bounds 
shown in the central panel of Fig. [2] are presently the 
strongest limit on the new-physics scale A. It is evident 
from the figure that including one-loop contributions im- 
proves the bound on inferred from the relic abun- 
dance by a factor of approximately 15. 

In summary, while collider bounds on DM-quark con- 
tact operators with Yukawa-like couplings are relatively 
weak when only tree-level processes are considered, much 
stronger bounds can be obtained by including heavy- 
quark loops in the analysis. In this letter, we used the 
recent CMS 5.0 fbT 1 and ATLAS 4.7 fbT 1 searches for 
jets with $ T to obtain the strongest collider limits on 
mass-dependent DM-quark scalar and pseudoscalar con- 
tact operators. Given that the LHC high-p^ experiments 
are rapidly accumulating luminosity, constraints on all 
contact operators involving DM particles with masses 
below the electroweak scale will improve significantly in 
the near future. The methods outlined here will be im- 
portant for further advances in constraining the parame- 
ter space of DM-quark interactions with mass-dependent 
couplings, and we are looking forward to their implemen- 
tation in future LHC analyses. 
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